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Abstract 

In this talk I discuss how a hrst order phase transition may proceed in rapidly expanding 
partonic matter produced in a relativistic heavy-ion collision. The resulting picture is that 
a strong collective flow of matter will lead to the fragmentation of a metastable phase into 
droplets. If the transition from quark-gluon plasma to hadron gas is of the first order, it 
will manifest itself by strong nonstatistical fluctuations in observable hadron distributions. 
I discuss shortly existing experimental data on the multiplicity fluctuations. 

1 Introduction 

A general goal of present and future experiments with heavy-ion beams is to study the 
properties of strongly interacting matter away from the nuclear ground state. The main 
interest is focussed on searching for possible phase transitions. Several phase transitions are 
predicted in different domains of temperature T and baryon density pb- As well known, 
strongly interacting matter has at least one multi-baryon bound state at Pb = Po ^ 0.16 
fm“^ corresponding to normal nuclei. It follows from the very existence of this bound state 
that there should be a first order phase transition of the liquid-gas type in normal nuclear 
matter at subsaturation densities, pb < Po, and low temperatures, T < 10 MeV. This phase 
transition manifests itself in a remarkable phenomenon known as nuclear multifragmentation. 

The situation at high T and nonzero baryon chemical potential pb {.Pb > 0) is not so 
clear, although everybody is sure that the deconflnement and chiral transitions should occur 
somewhere. Reliable lattice calculations exist only for /x = 0 i.e. ps = 0 where they predict a 
smooth deconflnement transition (crossover) at T 170 MeV. As model calculations show, 
the phase diagram in the (T, ps) plane may contain a first order transition line (below 
called the critical line) which ends at a (tri)critical point jHElE]- Possible signatures of 
this point in heavy-ion collisions are discussed in ref. |1]. However, it is unclear at present 
whether critical fluctuations associated with the second order phase transition can develop 
in a rapidly expanding system produced in a relativistic heavy-ion collison because of the 
critical slowing down effect [5]. In my opinion, more promissing strategy would be to search 
for a first order phase transition which should have much more spectacular manifestations 
jH] discussed below. It is interesting to note that, under certain non-equilibrium conditions, 
a first order transition is also predicted for symmetric quark-antiquark matter with zero net 
baryon density [7j. 

A striking feature of central heavy-ion collisions at high energies, confirmed in many 
experiments (see e.g. mi), is a very strong collective expansion of matter at later stages 
of the reaction. This process looks like an explosion with the matter flow velocities com¬ 
parable with the speed of light. The applicability of equilibrium concepts for describing 


phase transitions under such conditions becomes questionable and one should expect strong 
non-equilibrium effects uni im ini. In this talk I demonstrate that non-equilibrium phase 
transitions in rapidly expanding matter can lead to interesting phenomena which, in a certain 
sense, are even easier to observe. 

2 Effective thermodynamic potential 

To make the discussion more concrete, in this talk I adopt a picture of the chiral phase 
transition predicted by the linear sigma-model with constituent quarks 0 [ini. Then the 
mean chiral held $ = {a, n) serves as an order parameter. The constituent quark mass is 
generated by interaction with the sigma held, m = ga, where 5 ^ is a corresponding coupling 
constant. The ehective thermodynamic potential T, /r) depends, besides <h, on tempera¬ 
ture T and quark chemical potential /i = /is/3. This model respects chiral symmetry which 
is spontaneously broken in the vacuum, where the sigma held has a nonzero expectation 
value, (cr) = f-^, (tt) = 0. It is important for our discussion below that the model predicts a 
phase diagram on the (T,/i) plane with a critical point at (T=100 Mev. /i=207 MeV) and 
a hrst order phase transition line at lower T and /i. A schematic behaviour of r2(T,/i; $) as 
a function of the order parameter held cr at tt = 0 is shown in Fig. 1. The minima of hi 
determine the stable or metastable states of matter under the condition of thermodynamical 
equilibrium, where the pressure is P = —Vlmin/V. The curves from bottom to top corre¬ 
spond to homogeneous matter at diherent quark chemical potentials and hxed temperature 
T = 0. The dash-dotted curve corresponds to the hrst order phase transition point (two 
equal minima separated by a potential barrier). Two dashed curves show the thermody¬ 
namic potential at upper and lower spinodal points, where one of the minima disappears. 
The range of thermodynamic parameters where two phases, one stable and one metastable, 
may exist simulteneously is constrained by these two curves. The critical point would cor¬ 
respond to the situation when two minima fuse and the barrier disappears. This situation 
is illustrated in Fig. 2 showing the thermodynamic potential at diherent temperatures and 
hxed chemical potential fic = 207 MeV. This model reveals a rather weak hrst order phase 
transition, although some other models mis predict a stronger transition. The discussion 
below is quite general. 

One can plot a family of curves for the (T, /i) values corresponding to an isentropic 
expansion of matter. Qualitatively the potential curves look similar to the ones depicted 
in Fig. 1 . Assume that at some early stage of the reaction the thermal equilibrium is 
established, and partonic matter is in a “high energy density” phase Q (lowest curve). This 
state corresponds to the absolute minimum of O with the order parameter close to zero, 
(T ~ 0, TT = 0, and chiral symmetry restored. Due to a very high internal pressure, Q matter 
will expand and cool down. At some stage a metastable minimum appears in D at a hnite 
value of a corresponding to a “low energy density” phase H (lower dashed curve), in which 
chiral symmetry is spontaneously broken. At some later time, the critical line in the (T, //) 
plane is crossed where the Q and H minima have equal depths, i.e. Ph = -Pq (dot-dashed 
curve). At later times the H phase becomes more favorable, but the two phases are still 
separated by a potential barrier. At certain stage the minimum corresponding to the Q 
phase dissapears (upper dashed curve). The dashed curves separate the regions in the phase 
diagram where one of the phases is unstable (spinodal points). If the expansion of the Q 
phase continues until the barrier vanishes, the system will hnd itself in an absolutely unstable 
state at a maximum of the thermodynamic potential. Therefore, it will freely roll down into 
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Figure 1: The thermodynamical potential 02 for the sigma model at T = 0 and quark 
chemical potentials (starting from the top);/i = [0, 225, 279, 306, 322, 345, 375] MeV. 



Figure 2; The thermodynamical potential 02 for the sigma model for the sigma model at 
/i hxed to 207 MeV and temperatures (starting from the top):T = [0,50,75,100,125,150] 
MeV. 







the lower energy state corresponding to the H phase. This situation is known as spinodal 
instability. As shown e.g. in ref. jHj, the characteristic time of the ’’rolling down” process 
is relatively short, of about 1 fm/c. 

As well known, a first order phase transition proceeds through the nucleation process. 
According to the standard theory of homogeneous nucleation us, supercritical bubbles of 
the H phase can appear only below the critical line, when Ph > Pq- Then the critical 
radius for growing bibbles is = 2'j/{Ph — Pq), where 7 is the interface energy per unit 
area (surface tention). In rapidly expanding matter the nucleation picture might be very 
different. As argued in ref. jB], the phase separation in this case can begin as early as the 
metastable H state appears in the thermodynamic potential, and a stable interface between 
the two phases may exist. An appreciable amount of nucleation bubbles and even empty 
cavities may be created already above the critical line. They are stabilized by the collective 
expansion of matter. 

The bubble formation and growth will also continue below the critical line. Previously 
formed bubbles will now grow faster due to increasing pressure difference, Ph - Pq > 0 , 
between the two phases. It is most likely that the conversion of Q matter on the bubble 
boundary is not fast enough to saturate the H phase. Therefore, a fast expansion may lead to 
a deeper cooling of the H phase inside the bubbles compared to the surrounding Q matter. 
Strictly speaking, such a system cannot be characterized by the unique temperature. At 
some stage the H bubbles will percolate, and the topology of the system will change to 
isolated regions of the Q phase (Q droplets) surrounded by the undersaturated vapor of the 
H phase. 

3 Fragmentation of a metastable phase 

The characteristic droplet size can be estimated by applying the energy balance consider¬ 
ation, proposed by Grady nnum in the study of dynamical fragmentation of fluids. The 
idea is that the fragmentation of expanding matter is a local process minimizing the sum of 
surface and kinetic (dilational) energies per fragment volume. As shown in ref. [23, this pre¬ 
scription works fairly well for the multifragmentation of expanding nuclear systems produced 
in intermediate-energy heavy-ion reactions, where the standard statistical approach fails. 

Let us consider hrst an isotropically expanding system where the collective velocity held 
follows locally the Hubble law 

v{r) = H ■ r , (1) 

where H is a Hubble constant. Since there is no preferred direction, the droplets of the Q 
phase will have a more or less spherical shape. The characteristic radius of the droplet, R, 
can be estimated as follows ng. The energy AE associated with a Q droplet embedded in 
a background of the H phase is represented as the sum of three terms, 

AE = Ejjulk + Ekin + E SUV (2) 

The bulk term is simply equal to A£V where A£ = £q — £h is the energy density difference 
between the two bulk phases and H oc is the volume of the droplet. The second term is 
the collective kinetic energy of the droplet expansion with respect to its center of mass, 

Ekin = ^l £v\r)dV = ^A£VH^R\ 


( 3 ) 


The last term in eq. m is the interface energy, Egur = 7 * 5 ', which is parametrized in terms 
of the effective surface tension 7 and the surface area S oc R^. Grady’s argument is that 
the redistribution of matter is a local process minimizing the energy per droplet volume, 
AE/V. Then the bulk contribution does not depend on R, and the minimization condition 
constitutes the balance between the collective kinetic energy and interface energy. This leads 
to the optimum radius of a droplet 
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It is worth noting that the collective kinetic energy acts here as an effective long-range 
potential similar to the Coulomb potential in nuclei. 

At ultrarelativistic collision energies associated with RHIC and LHC experiments, the 
expansion of partonic matter will be very anisotropic with its strongest component along 
the beam direction ^Hj. Clear indications of such an anisotropy are seen already at SPS 
energies (see jH]). It is natural to think that in this case the inhomogeneities associated with 
the phase transition will rearrange into pancake-like slabs of Q matter embedded in a dilute 
H phase. The characteristic width of the slab, 2L, can be estimated in a similar way and 
the resulting expression for L* differs from Eq. 0 only by a geometrical factor (3 instead 
of 5 in parentheses). Generally, the faster is the expansion, the smaller are the fractures. Of 
course, at a later time the Q droplets will further fragment in the transverse direction due 
to the standard nucleation process. 

As Eq. 0 indicates, the droplet size depends strongly on H. When expansion is slow 
(small H) the droplets are big. Ultimately, the process may look like a hssion of a cloud of 
plasma. But fast expansion should lead to very small droplets. This state of matter is very 
far from thermodynamical equilibrium, particularly because the H phase is very dilute. One 
can say that the metastable Q matter is torn apart by a mechanical strain associated with 
the collective expansion. This has a direct analogy with the fragmentation of pressurized 
fluids leaving nozzles. In a similar way, splashed water makes droplets which have nothing 
to do with the liquid-gas phase transition. 

The driving force for expansion is the pressure gradient, VP == c^V£^, which depends 
crucially on the sound velocity in the matter, Cg. Here we are interested in the expansion 
rate of the partonic phase which is not directly observable. In the vicinity of the phase 
transition, one may expect a “soft point ” [131201 where the sound velocity is smallest and 
the ability of matter to generate the collective expansion is minimal. If the initial state of 
the Q phase is close to this point, its subsequent expansion will be slow. Accordingly, the 
droplets produced in this case will be big. When moving away from the soft point, one 
would see smaller and smaller droplets. For numerical estimates we choose two values of 
the Hubble constant: P“^=20 fm/c to represent the slow expansion from the soft point and 
H~^=6 fm/c for the fast expansion. 

One should also specify two other parameters, 7 and AS. The surface tension 7 is a sub¬ 
ject of debate at present. Lattice simulations indicate that at the critical point it could be as 
low as a few MeV/fm^. However, for our non-equilibrium scenario, more appropriate values 
are closer to 10-20 MeV/fm^ which follow from effective chiral models. As a compromise, 
the value 7 = 10 MeV/fm^ is used below. Bearing in mind that nucleons and heavy mesons 
are the smallest droplets of the Q phase, one can take AS = 0.5 GeV/fm^, i.e. the energy 
density inside the nucleon. Then one gets i?*=3.4 fm for H~^=20 fm/c and R*=1.5 fm for 
H~^=6 fm/c. As follows from eq. (|3)), for a spherical droplet V oc 1/AS, and in the hrst 



approximation its mass, 


(5) 


M* ^ ASV 
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is independent of AS. For two values of R* given above the mass is ~100 GeV and ~10 
GeV, respectively. The pancake-like droplets could be heavier due to their larger transverse 
size. Using the minimum information principle one can show II71I2II that the distribution 
of droplets should follow an exponential law, exp Thus, about 2/3 of droplets have 

masses smaller than M*, but with 1% probability one can find droplets as heavy as 5M*. 


4 Observable manifestations of quark droplets 


After separation, the droplets recede from each other according to the Hubble law, like galax¬ 
ies in expanding Universe. Therefore, their c.m. rapidities are in one-to-one correspondence 
with their spatial positions. One may expect that they are distributed more or less uniformly 
between the target and the projectile rapidities. On this late stage it is unlikely that the 
thermodynamical equilibrium is re-established between the Q and H phases or within the 
H phase alone. If this were to happen, the hnal H phase would be uniform, and thus there 
would be no traces of the droplet phase in the hnal state. 

The hnal fate of individual droplets depends on their sizes and on details of the equation 
of state. Due to the additional Laplace pressure, 27 /i?, the residual expansion of individual 
droplets will slow down. The smaller droplets may even reverse their expansion and cooling 
to shrinking and reheating. Then, the conversion of Q matter into H phase may proceed 
through the formation of the imploding dehagration front 1201E2]. Bigger droplets may 
expand further until they enter the region of spinodal instability At this stage the diherence 
between 1-st and 2-nd order phase transitions or a crossover is insignihcant. Since the 
characteristic “rolling down” time is rather short, ~ 1 fm/c, the Q droplets will be rapidly 
converted into the non-equilibrium H phase. In refs. uni ini the evolution of individual 
droplets was studied numerically within a hydrodynamical approach including dynamical 
chiral helds (Ghiral Fluid Dynamics). It has been demonstrated that the energy released at 
the spinodal decomposition can be transferred directly into the collective oscillations of the 
{a, it) helds which give rise to the soft pion radiation. One can also expect the formation of 
Disoriented Ghiral Gondensates (DGG) in the voids between the Q droplets. 

An interesting possibility arises if the metastable Q phase has a point of zero pressure. 
In particular, this is the case for the MIT bag model equation of state at temperatures only 
slightly below In this case the droplets might be in mechanical equilibrium with 

the surrounding vacuum {Ph ~0), like atomic nuclei or water droplets. The equilibrium 
condition is 


Pq = 
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( 6 ) 


where vq = 12 is the degeneracy factor for massless u and d quarks (the gluon contribution 
is omitted here), and H is a bag constant. The evolution is then governed by the evaporation 
of hadrons from the surface (see also the discussion in Ref. jUj). One can speculate about 
all kinds of exotic objects, like e.g. strangelets, glueballs, formed in this way. The possibility 
of forming ’’vacuum bubbles’, i.e. regions with depleted quark and gluon condensates, was 
discussed in ref. HD]. All these interesting possibilities deserve further study and numerical 
simulations. 
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Figure 3: Schematic view of the momentum space distribution of secondary hadrons pro¬ 
duced from an ensemble of droplets. Each droplet emtts hadrons (mostly pions) within a 
rapidity interval 6y ^ 1 and azimuthal angle spreading of 6(j) ~ 1. 

After separation the QGP droplets recede from each other according to the global expan¬ 
sion, predominantly along the beam direction. Hence their center-of-mass rapidities yi are 
in one-to-one correspondence with their spatial positions. Presumably yi will be distributed 
more or less evenly between the target and projectile rapidities. Since rescatterings in the 
dilute H phase are rare, most hadrons produced from individual droplets will go directly into 
detectors. This may explain why freeze-out parameters extracted from the hadronic yields 
are always very close to the phase transition boundary 123 . 

In the droplet phase the mean number of produced hadrons in a given rapidity interval 
is 

No 

(N) = '^fk = {n){ND) , (7) 

i 

where nl is the mean multiplicity of hadrons emitted from a droplet i, (n) is the average 
multiplicity per droplet and {Nd) is the mean number of droplets produced in this interval. 
If droplets do not overlap in the rapidity space, each droplet will give a bump in the hadron 
rapidity distribution around its center-of-mass rapidity yi ^3]. In case of a Boltzmann 
spectrum the width of the bump will be dy ~ ^JT/m, where T is the droplet temperature 
and m is the particle mass. At T ~ 100 MeV this gives dy ~ 0.8 for pions and 5y ~ 0.3 for 
nucleons. These spectra might be slightly modihed by the residual expansion of droplets. 
Due to the radial expansion of the hreball the droplets should also be well separated in the 
azimuthal angle. The characteristic angular spreading of pions produced by an individual 
droplet is determined by the ratio of the thermal momentum of emitted pions to their 
mean transverse momentum, 5(j) ~ 3T/(p^) ~ 1. The resulting phase-space distribution of 
hadrons in a single event will be a superposition of contributions from different Q droplets 
superimposed on a more or less uniform background from the H phase. Such a distribution 
is shown schemmatically in Fig. 3. It is obvious that such inhomogeneities (clusterization) 
in the momentum space will reveal strong non-statistical fluctuations. The fluctuations will 
be more pronounced if primordial droplets are big, as expected in the vicinity of the soft 















point. If droplets as heavy as 100 GeV are formed, each of them will emit up to ~200 pions 
within a narrow rapidity and angular intervals, hr] ~ 1, 50 ~ 1. If only a few droplets are 
produced in average per unit rapidity. No ^ 1, they will be easily resolved and analyzed. 
On the other hand, the fluctuations will be suppressed by factor ^/ND if many small droplets 
shine in the same rapidity interval. 


5 Anomalous multiplicity fluctuations 


For our discussion below we consider a more general case of the droplet mass distribution 
when masses follow a gamma-distribution 

Wk{M) = exp {-bM) , (8) 

T{k) 

which is normalized for 0 < M < oo. The mean mass and its standard deviation are 
expressed through the parameters k and b as 


, , k ^/k 

< M >= - , aM = -r- 
b b 
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^/k 


(9) 


These expressions show that quantity 1/y/k gives the relative scale of fluctuations of M 
around < M >. It should be stressed that the gamma-distribution (jS)) drops at large M less 
rapidly than a corresponding gaussian distribution. 

One can easily calculate the combined multiplicity distribution produced by the ensemble 
of many droplets. Let us assume that the normalized mass distribution of droplets is Wk{M) 
and that each droplet emits hadrons according to the Poisson law, Pn{n), with the mean 
multiplicity proportional to the droplet mass, n = M/{Et,) (for pions {Et,) ^ 3T 0.5 
GeV). Then the combined distribution is given by the convolution of the two, 

f-M 

Pk{N) = / dMwk{M)p^{N). (10) 

Jo 

For the gamma-distribution (jH)) one can perform explicit analytical calculations. It is re¬ 
markable that the resulting distribution is a famous Negative Binomial Distribution (NBD) 
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( 11 ) 


In a limiting case of the exponential mass distribution (k=l) the combined distribution is 
simply given by 
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( 12 ) 


where (N) = (n) = a{M) is the mean total multiplicity. 

It is convenient to characterize the fluctuations by the scaled variance 
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( 13 ) 








Its important property is that a;Ar = 1 for the Poisson distribution, and therefore any devi¬ 
ation from unity will signal a non-statistical emission mechanism. For the NBD, eq. (??), 
one easily hnds ujn = 1 + {N)/k. As shown in ref. [21], for an ensemble of emitting sources 
(droplets) a ;at can be expressed in a simple form, cuat = a;„ -f {n)ujD, where Un is an av¬ 
erage multiplicity fluctuation in a single droplet, ud is the fluctuation in the droplet size 
distribution and (n) is the mean multiplicity from a single droplet. Since a;„ and ojd are 
typically of order of unity, the fluctuations from the multi-droplet emission are enhanced by 
the factor (n). According to the picture of a hrst order phase transition advocated above, 
this enhancement factor could be as large as 10^. It is clear that the nontrivial structure 
of the hadronic spectra will be washed out to a great extent when averaging over many 
events. Therefore, more sophisticated methods of the event sample analysis should be ap¬ 
plied. As demonstrated below, the simplest one is to search for non-statistical fluctuations 
in the hadron multiplicity distributions measured in a varied rapidity bin. 



Figure 4: Event-by-event multiplicity distributions in dr] intervals 0.1, 0.2, 0.3, 0.5, 1.0 
measured by E-802 Collaboration for central ^®0-|-Cu collisions at lab energy 14.6 AGeV 
P7j . The data for each interval are plotted as a function of n/{n) and scaled by (n), the 
mean multiplicity in the interval. Each successive distribytion has been normalized by the 
factor indicated in the hgure. The shape evolvs from almost Gaussian {6t] = 1.0) to nearly 
exponential (577=0.1). 

The event-by-event multiplicity fluctuations observed so far in heavy-ion experiments 
do not reveal any anomalous enhancement. Fig. 4 shows an example of the multiplicity 
distributions in varying pseudorapidity intervals measured for 0-l-Cu central collisions at 
AGS j2Zl- These distributions can very well be htted by the NBD where parameters k{ST]) 
and {N{6r])) follow a linear relationship with nonzero intercept k{0). The data show an 
increase in scaled variances of about 15-i-30% over 1. Apparantly such moderate deviations 
can be explained by ordinary reasons, not related to a phase transition. An interseting 
observation has been made by the NA49 collaboration j2Hl; which has found non-monotonic 
behaviour of un as a function of the projectile participant number Np. But in this case too 
the actual values of a ;at are only of about 2, and most likely they can be explained by the 
fluctuations in the number of target participants 


















It should be noted that the NBD hts were hrst used to describe the multiplicity distri¬ 
butions in high-energy pp and pp collions (see e.g. ref. ffl). They are consistent with the so 
called KNO scaling 


6 Conclusions 

• A hrst order phase transition in rapidly expanding matter should proceed through 
the nonequilibrium stage when a metastable phase splits into droplets whoose size is 
inversly proposrtional to the expansion rate. The primordial droplets should be biggest 
in the vicinity of a soft point when the expansion is slowest. 

• Hadron emission from droplets of the quark-gluon plasma should lead to large nonsta- 
tistical huctuations in their rapidity and azimuthal spectra, as well as in multiplicity 
distributions in a given rapidity window. The hadron abundances may rehect directly 
the chemical composition in the plasma phase. 

• To identify the phase transition threshold the measurements should be done at different 
collision energies. The predicted dependence on the expansion rate and the reaction 
geometry can be checked in collisions with different ion masses and impact parameters. 

• If the hrst order deconhnement/chiral phase transition is only possible at hnite baryon 
densities, one should try to identify it by searching for the anomalous huctuations in 
the regions of phase space characterized by a large baryon chemical potential. These 
could be the nuclear fragmentation regions in collisons with very high energies (high- 
energy SPS, RHIC, LHC) or the central rapidity region (AGS, low-energy SPS, future 
GSI facility FAIR). 

The author is grateful to L.M. Satarov, M.I. Gorenstein and M. Gazdzicki for many fruit¬ 
ful discussions. This work was supported partly by the RFBR grant 05-02-04013 (Russia). 
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